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Abstract 

Gravitational backgrounds md+2 dimensions have been proposed as holo- 
graphic duals to Lifshitz-like theories describing critical phenomena in d + 1 
dimensions with critical exponent z > 1. We numerically explore a dilaton- 
Einstein-Maxwell model admitting such backgrounds as solutions. Such back- 
grounds are characterized by a temperature T and chemical potential /j,, and 
we find how to embed these solutions into AdS for a range of values of z and 
d. We find no thermal instability going from the T <^ /ito the T 2> regimes, 
regardless of the dimension, and find that the solutions smoothly interpolate 
between the Lifshitz-like behaviour and the relativistic AdS-like behaviour. 
We exploit some conserved quantities to find a relationship between the en- 
ergy density £, entropy density s, and number density n, £ = {Ts + fin), 
as is required by the isometries of AdS(i+2- Finally, in the T <^ jjL regime the 
entropy density is found to satisfy a power law s oc cT'^/^^u^^"^)^/^, and we nu- 
merically explore the dependence of the constant c, a measure of the number 
of degrees of freedom, on d and z. 



1 Introduction 



By now, holography is considered a general tool for studying many kinds of physical 
systems. In extremely symmetric situations, e.g. the AdSsX background where 
AdS/CFT was born [l], the symmetries allow many qualitative and quantitative 
checks of the duality. A vast number of setups for superconformal field theories 
have been studied |2), again taking advantage of the symmetries. Less symmetric 
situations exist [s], but these still are designed to address questions relevant for 
particle physics, and so preserve Poincare symmetry by construction. 

Much has been learned from the 3 + 1 dimensional cases described above. It is 
known that finite temperature in the field theory often corresponds to the presence 
of a horizon in the gravitational dual. Further, U{1) gauge symmetries in the bulk 
correspond to conserved number operators in the dual field theory. Therefore, to 
study a field theory at finite temperature and chemical potential in a holographic 
setup, the focus is on charged objects with horizons in the bulk: charged black 
holes H]. 

Recently these techniques have been brought to bear on other types of systems. 
Most notably, condensed matter systems have been a topic of intense research us- 
ing holographic techniques [5]-[8] (see |9 11 for an introduction). This seems quite 



natural, given that theoretical condensed matter and particle physics share quantum 
field theory (QFT) as their primary language. The intuition gained from the 3+1 
dimensional setups suggest that holography should also be useful to study condensed 
matter systems at strong coupling, further implying that certain condensed matter 
systems may have regimes of "stringy" behaviour. 

In particular, much effort has gone into describing quantum critical behaviour of 
these theories. As mentioned above, identifying symmetries is of primary importance. 
Quantum critical systems typically exhibit a scaling symmetry 

t A^t, Xi Xxi (1) 

similar to the scaling invariance of pure AdS (2; = 1) in the Poincare patch. From a 
holographic standpoint, this suggests the form of the spacetime metric 

ds^ = (^r^'dt^ + r^dx'dx^ij + , z G {1, 2...d} (2) 

where the above scaling is realized as an isometry of the metric along with r — t- A^^r. 
The above metric has no symmetry that mixes time and space, although there are 



also models that contain a Galilean (or other non-relativistic) symmetry |5lp, 12 



many of which admit (or were constructed using) supergravity embeddings. It is 
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interesting to note that many of the metrics used for non-relativistic physics are 



known to be coset spaces 13 



While it may at first appear as though g is a simple generalization of AdS, 
there are some significant differences |14j. Here we will not address such concerns, 
but it should be noted that the above spacetime has been embedded into lOD and 



IID supergravity theories 15[j21 and it would be interesting to see how, or if, these 

Further these models 



theories avoid the problems discussed in 14 



15 19 21 



may 



shed light onto the weak coupling degrees of freedom, and so show what types of 
field theories these gravitational systems describe. 

Here we will be concerned with phenomenological models that admit the metric 
([2]) as a solution. Two such models with d = 2 have Lagrangians given by 22 



S' 
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2A - -g^ A 
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(3) 



(we call this model S' for short) and 23 



^ = TTT?^ / ^'^v^ (i? - 2A - 2(V0)2 - e'-^g') , 
lo7rG4 J 



(4) 



(we call this model S for short) where in either case g = dA is a two form field 
strength. There are also examples of actions with corrections that admit solutions 
of the form ^ [24]. 

There are several differences between the above models. First, the Lifshitz solu- 
tion to model S' is truly invariant under the Lifshitz rescaling symmetry, while the 
solution with metric ^ of model S has a logarithmically running dilaton, and so we 
call this a "Lifshitz-like" geometry. Note also that S has a [/(I) gauge symmetry 
while S' does not, and so the field theory dual to S has a conserved particle number 
with chemical potential /i. Also of note is that model S admits an exact black 
brane solution 23 that asymptotes to the metric (|2]), which can also be generalized 
to higher dimensions 25 . Finding black brane solutions for action 5" has proven 
more difficult, and one often needs to resort to numeric methods [26-29 , but not 
always 27, 30] (an analogous analytic statement for an R^ extension may be found 
31 , and one should also see l32] where a certain extension to this model admits an 



analytic black hole). For extensions and variations to the basic model S", see 31 ,33 



and for the holographically renormalized action, see 34 
Similar actions to S can be found in 



35 -46 , which either contain different matter 



fields (probe, or back- reacted) , different couplings, or some combination. However, 
here we will study the action S, also studied in |23l|29l|47|48|, generalized to arbitrary 
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dimension as in 



25 



S = \ [ d^+^x^ (R-2A- 2(V</))' - e^'^'^g^) . (5) 
167rGrf+2 J 

We will consider "UV completing" by embedding black branes into asymptotically 
AdS space, generalizing to generic d the discussion of |48j. We will explore the 
thermodynamics of these backgrounds analytically and numerically. 

As investigated in our previous work |48|, there are two scales in the problem, 
T and /z, and there are two regimes, T <^ /i (Lifshitz-like) and T ^ ^ (AdS-like). 
Therefore, one might anticipate a thermal instability during the transition between 
the two regimes. However, in the d = 2 case, we found that there was no such 
instability, and so no discontinuous phase transition. We thought that this could be 
related to a Coleman— Mermin— Wagner theorem applied to the thermal vacua of a 
2 + 1 dimensional theory. However, here we find that there is no discontinuous phase 
transition, and the model smoothly interpolates between the Lifshitz-like behaviour 
and the AdS-like behaviour, regardless of d. 

Finally, in the Lifshitz-like regime, we expect a relation of the form 

4Gd+2S = c{z,d){LfiY(^^^ ' , (6) 

where s is the entropy density, T is the temperature, d is the number of spatial 
dimensions, and z is the critical exponent. The coefficient c(z, rf) is in some sense a 
measure of the number of degrees of freedom of the system, and we find its behaviour 
function of d and z numerically. 
The rest of the paper is organized as follows. In Section |2] we reduce the model 
to an effective radial model for arbitrary dimension rf, and find that the equations of 
motion reduce to a set of 4 first order ordinary differential equations. We comment 
on the proper normalization of the charge density q. We compute the perturbative 
expansion at the horizon and at the AdS asymptotic, and use this to show the relation 

£._l_,r. + ,„) (7) 



which is necessary by the scaling symmetry of AdS 48 . 

Finally, in Section |3j we turn to numeric results. We set up the numerical inte- 
gration, and show how to efficiently parameterize the different regimes of T j fi with 
one piece of data specified at the horizon. We then numerically integrate, and show 
that for d = 2,3,4- ■ -9 that there is no discontinuous phase transition: the solu- 
tions smoothly and monotonically interpolate between the two regimes T/fi <^ 1 
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and T//i ^ 1. We then find c{z, d) for a range of z and and discuss some of its 
qualitative behaviour. 



2 Analysis of the model 
2.1 Reduction 

We wish to consider dilatonic black brane solutions to the equations of motion fol- 
lowing from the action 

S = \- — I d'^+^x^/^ (R-2A- 2(V0)2 - e^'"''g^) . (8) 

16TTGd+2 J 

We reduce the above action on the following Ansatz 



ds^ = -e^^('-)dt' + e^^^'-) (dx'dx^,,) + e^^^'^dr 

= 0(r), A = e^^''^dt (9) 

where i G {1,2, ■ ■ ■ d} and Q = dA. We will get only ordinary differential equations 
in what follows, and so we define d = We may reduce the action to a one 
dimensional action Q 

S = ^^J^^ ^ [ 2dtl[dx, I drLiD (10) 



:+2 



i=l 

where the one dimensional Lagrangian is given by 

Note that when d = 2 this Lagrangian agrees with our earlier work. It can be verified 
that all equations of motion associated with the action ([s]) are reproduced by (11 ), as 
long as one uses the equation of motion for C. Here, C acts as a Lagrange multiplier, 
imposing the "zero energy" condition. Further, we note that C allows for generic 
r diffeomorphisms. We will refer to changing r coordinate as "coordinate gauge" 
transformations, to differentiate from the U{1) gauge transformations associated with 
A. 

^We keep track of normalization for later use. 
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There are many conserved quantities associated with the above action. First, 
there is the conserved quantity associated with the shift symmetry {A, B,C,(f),G) — >■ 
{A + dSi, B — Si,C,(l),G + dSi). This can be understood as a rescahng of the time 
coordinate, and the Xi coordinates that leaves dtYlidxi invariant, and this then 
descends to the ID action as Noether symmetry. Further, we note that there is the 
conserved quantity associated with e*^ — > e'^ + const which is exactly the global part 
of the gauge symmetry associated with A. Further, we note that (^4, B,C,(f),G) — >■ 
{A, B,C,(I) + S2, G — aS2) is also a symmetry. One can view this as saying that the 
dilaton here is a space dependent gauge coupling, and we may absorb a normalization 
of the gauge coupling into the definition of A. In addition, we have the Hamiltonian 
constraint as always. Now we count 4 dynamical fields, and 4 conserved quantities. 
This imphes that the system is completely equivalent to a set of first order differential 
equations 

_^A+dB-c^Q^y + e^+'^^+^A = 0, (12) 

^^A+dB-Cg^ _ ^^A+dB-CQQ _ ^^^-A+dB-C+2G+2a^QQ ^ (^3) 
^-A+dB-C+G+2a^gQ = Q. (15) 

We form the combination Dq — uVq + dVo, Pq — \Vo and find 

de^+'^^-CdAdB + d{d-l) ^A+dB-C^QQY ^ ^^A+dB-C+2G+2.^^QQY 

_^A+dB-c^Q^Y + e^+'^^+^A = 0, (16) 
d{2d(l) + a{dA-dB))e^+'^^~^ = Df,, (17) 
^A+dB-CQ^ + ae-'^+'^^-^+2«+^"^aG = Po, (18) 

^.A+dB-C+G+2a4.QQ ^ Q (^g) 

There are two known one parameter families of solutions. First, there is the AdS 
black brane 



A{r) = In I Lry^l - (J^Y^^ j , B{r) = In(Lr), 
C(r) = In 



1 - m'^' 

A = Qbdt. (20) 
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where A = —d{d — 1)/(2L^), and gt and (pb are arbitrary constants. For this back- 
ground, the conserved quantities are Q = 0, Dq = d{d + l)L'^ar'l'^^ /2 and Pq = 
0. Further, the solution has T = LT = rh{d + 1) / (Att) by reading off the pe- 
riodicity of imaginary time at the horizon. Further, using the area law for the 
entropy density, we find that s = rf/{4Gd+2), or in thermodynamic terms s = 
{47ryT'^L^/[{d + lY4Gd+2] ■ This shows that the number of de grees of freedom in 
the field theory is linear in L'^/[AGd+2)- 



There is also the Lifshitz black brane given by 23 , 25 



A{r) = In [ LLttLr'Jl- (^j"^'^ ) , B{r) = ln(Lr 



C(r) = In — ^ , 2a0(r) = In (r-^'^*) 



r 



1 _ (^!iLy+'^ 



G(r) = In ^ '- I . (21) 

The constants are given by 

J_z + d){z + d-ll 2{Q/L'^-')\z + d-l) 

^^"^ d{d+l) ' d{d + l){z-l) ' ^^^^ 

ai is arbitrary (and may be removed by time rescaling), and a = >y2d/{z — 1). The 
conserved quantities are Pq = 0, Dq = dj \l(z — V)\r^y^'^aLlj'^d[z^ d) and Q (which 
is arbitrary). Given that i?(r) is normalized in the same way as the AdS black brane, 
we read s = rjJ/[4Grf+2]- Finally, we see that T = r\fli^[z -|- (i)/(47r). While this is a 
dimension free measure of temperature, it is not clear what units to use. This is to 
be expected for a non-relativistic theory because energies and length scales are not 
interchangeable. However, this ambiguity is only a constant, and so s oc T'^l\ When 
we consider embedding these black branes into AdS space the ambiguity is removed. 



2.2 Perturbation theory at the horizon 

We first expand the functions as "correction functions" that should go to constant 
values for an AdS background. We therefore define 

yl(r) = In(Lr) + y4i(r), 5(r) = In(Lr), 

C(r) =ln f-) +Ci(r), G(r) = ln(L) + Gi(r). (23) 
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We further define dimensionless conserved quantities 

Q^L^-^Q, Do^L'^Do, Po^L^Pq^O (24) 

where we remind the reader that we will only use the Pq — gauge. We have 4 first 
order differential equations, but we only have 3 dynamical functions. Therefore, we 
may eliminate one function using an algebraic expression, and we do so for the field 

e-^'^^ - ^ I dr^'^+^aid + 1) (e^^+^i - e^^-^^ 

-2r^+^bo - Adr'^+^e^'aQ + 2ae^'-^' {aQe^'f^ . (25) 

Using this, the other differential equations may be written 

e^' (Do + 2daQe^A 
S-'' ^=0. (215) 

^ ^ = 0. (27) 

^ 2d+2 / 

9e«i - ^^^^^ {{d+1) (e^i+^i - e^-^i) (28) 



2aQr 



-2r'^+^ (^Do + 2daQe^'^ + 2ae^^-^^a^Q 



These functions may be expanded around a regular horizon as 

e^i = oo (ir - r^)^ + ai(r - rft)i H ) , 

„Ci _ Co 



- + ci(r-rft)2 ••• , 
= ao (^o(r - ru) + ^i(r - rft)^ + •••). (29) 
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One may read that the dilaton goes to a constant. With Pq = 0, we find that the 
equations of motion require 

2coDo _ a'^d{d + ifc^ - 2drh[c? - 2)(rf + \)cl + [o?d - 4 - 6rf) 

_ Co {Zc^did + \fcl - 2rfr^(2 + -ia^^d + l)c2 + rl (Sa^rf + 4 + 6rf)) 
Co — Co, Ci — — ^ , 

d{{d^\)cl-r^)ri 
91 = ( rt-'a\d + l)^c^ - ri-'a\d + 1)^0^ 



8(5 cor^ 



+ 2){d+l)cl + ri^'{a' + 2) . 



This will provide initial conditions for the equations of motion when we numerically 
integrate. From this expansion we can read the temperature and entropy density 

^ rla^ ^ Do ^ 

AnLco adrf2^V AGd+2 

Further, we find using the above type of expansion that the Lifshitz black brane satis- 
fies Co = \/ + 2^Jrh/ {d + 1) while the AdS black brane satisfies Co = ^/rh/ {d + 1). 
We expect that these two limits bound the physical values of Cq. 



2.3 Perturbation theory at r = oo: AdS asymptotics 



We want to expand the solution about infinity and require that it asymptotes to the 
pure AdS solution. For this, we use the expansion: 



A{r 
B{r- 



G{r 



In(Lr) + Ai{r) 
In(Lr), 



In 



+ Ci(r), 



(32) 



ln(<l>b) + 0i(r), 
In(^b) + Gi(r), 



where now Ai{r), Ci(r), G'i(r), and 0i(r) are understood to be perturhativ eiunciions. 
Inserting equations (32) in the equations of motion (16)-(19) and performing the 
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integrations gives (to leading order) 

Do + 2daQgb - 2dPo 



Ci{r) 
Gi{r) 

Mr) 



Do + 2daQgb - 2dPo 
d{d + l)aL'^r'i+^ 

Q 

(d- l)^5$2"L'^-2r'^-i^ 
aQgb - Po 



+ 



ci(ci- l)<l>2"L2{'i-i)r2'i' 

2Q2 

(d- l)(rf+ l)<l>fL2('i-i)r2««' 



{d + l)L'^r'^+i 2d{d - l)$2«L2('^-i)r 



2d' 



(33) 
(34) 
(35) 
(36) 



The constants of integration of equations ([TS]) and (19) (which would lead to con- 
stant modes in (36) and (35) ) are absorbed in the boundary values of $6 and gb , 
respectively. The constant of integration of equation (17) is removed by rescaling 
the time coordinate. 

We next want to evaluate the energy density. We follow closely our previous 
discussion for the d = 2 case 28 , and generahze it here to arbitrary d. We use the 
background subtraction technique in 49 , where the total energy density is given by 

1 

"8^ 



{N, {'K Ko) - Ntp,.r) 



(37) 



where Nf is the lapse function, Nj^ is the shift vector, '^K is the extrinsic curvature 
of the d-dimensional spatial boundary slice inside the constant-t slice, '^Ko is the 
(i-dimensional extrinsic curvature of the reference background, p^i, is the momentum 
conjugate to the time derivative of the metric on the constact-t slice, and f is the 
spatial unit vector normal to the constant-r surface. The reference background we 
take is pure AdS in the Poincare patch. 



For our metric, Nt 
£ = lim - 



0, f = e , and the energy density is 



2B\ Jd~2)B 



-Cri 







(38) 



For the pure AdS reference, we have e'^^'^f = ^. The total energy density then reads 

£ = lim -C^ , (39) 



lim -—^ -Ci , 



where we have dropped terms that go to zero in the large r limit. Therefore, we read 
the energy density to be 

1 2 {Do + 2daQgb) 



£ 



167rG'd+2 (rf+ 1) aL'i-^^ 



(40) 
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Now we need to normalize the charge density Q. For this, we follow the discussion 
of 50 , where they use the action to normalize Q in such a way that it is conjugate 
to the fixed field value at infinity. In our context, this reads 

Sg2 = Y^~^Y^l9eom, (41) 

where qgeom is a normalized charge density and figeom is the potential difference 
between the horizon and the boundary. For us, we are using a gauge where Pq = 
and so the potential at the horizon is 0, and so = ftgeom for this gauge choice. We 
calculate Sg2 and find 

^--16^b^^«' '''' 
and so we identify the charge Q as 

Q = qgeom- (43) 



Inserting the scaled quantities (24), and using the relations (31) for the temperature 
and entropy density along with the normalization of Q, we find that 

£ = ^ (T. + /xn) , (44) 

where we have changed to field theory quantities figeom = l^L'^, qgeom = n/L"^. This 
is the expected relation for a conformal field theory. 



2.4 Other considerations, and setup for numerics 

Given the discussion of perturbation theory around the horizon, we find that the 
constant cq has the following window of allowed values 

where the lower limit gives the pure AdS black brane, and the upper limit gives the 
pure Lifshitz black brane condition. Using this, we may efficiently parameterize all 
black branes in terms of this one quantity, which we now explain. 

First, we may rescale r by a constant: this will affect the location of the horizon, 
and effectively allows us to fix the horizon to be at r/^ = 1. After doing so, we may 
use time rescaling to set Dq to be any value we wish, and further we may use the 
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global symmetry associated with shifting the dilaton to move Q to be any value we 
wish. Therefore, we may fix all of the horizon data except for cq using symmetries. 
Then, cq parameterizes the different black branes, and the allowed region of cq is 
given by the above range (45). However, now the asymptotic values of the fields 
are non-canonical. This means that, given a Cq, we take these asymptotic values of 
the fields e^,e'^, to be the output values of the numeric integration started at the 
horizon. These output values encode how to use the time rescaling and the global 
symmetry to bring them to their canonical values. Hence, these parameterize a given 
T and /i for a fixed Vh = 1 black brane. Near the limiting values of the range (45) we 
expect the ratio T//i to go to zero (Lifshitz-like regime) or to infinity (AdS regime), 
and so all possible ratios of T/ fi are explored. 

To find the generic black brane with temperature T and chemical potential fi, 
we would simply find the appropriate black brane at r/j = 1 with the same value of 
T/ii. We would then use the scaling symmetry of AdS to adjust, say T, to its correct 
value. The resulting black brane has the specified value of T as its temperature, 
as its chemical potential, and the position of the horizon Vh will be determined by 
the rescaling. These will all be unique as long as the ratio T/fi occurs only once for 
the r/i = 1 reference black brane. The crucial question is then whether the graph 
of T/fi vs. Co is monotonic for the r/j = 1 black branes. If it is, there is always a 
unique black brane given a particular value of T and /i, and so there are no possible 
discontinuous phase transitions. Indeed, we find this is the case. For the sake of 
clarity, instead of graphing the quantity T/fi as a function of the somewhat esoteric 
quantity cq, we simply graph s vs. T for fixed /i. We will see that it is monotonic, 
and interpolates between the Lifshitz-like scaling s oc T'^/^ and the AdS behaviour 
s oc T'^. Accordingly, there is no discontinuous phase transition. 

Next, we recall a subtlety that emerged in 48 associated with the definition of 
the chemical potential /i, and we recall the resolution here. First, the physical value 
of 6*^1 is not strictly determined. We may use the global symmetry {A, B,C,(j),G) — ?■ 
{A, B,C,(f) + 62, G — aS2) to rescale this value to any value we wish. Because this is 
a global symmetry that does not involve the metric, the stress energy tensor (and 
therefore the geometry), is not determined by this number. In fact, only global 
symmetry invariants can determine anything in the geometry (this was also noticed 
in 47 for extremal solutions). We wish to consider /i as a scale in the theory, i.e. it 



is the scale at which new particles can be added/excited, and so we expect this to 
correspond to some scale in AdS, and so must be a global symmetry invariant. 

One way to fix this ambiguity is to note that the gauge coupling at infinity is the 
asymptotic value of e"^""^ (which also transforms under the global symmetry). Hence, 
we would like the gauge kinetic term to go to a canonical value. This can be arranged 
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by using the global symmetry to scale e'^°"^\r=oo = 1- From the beginning we should 
have expected such a statement: only fixing both the asymptotic value of e°"^ and 
the asymptotic value of e*^^ will determine the geometry. We had traded this for e*^^ 
and Q in the above discussions. Therefore, one can read the output value of /i from 
the numeric integration as 



(46) 



The factor of e^'^ in the denominator accounts for the fact that e*^^ transforms under 
the time rescaling used to set the asymptotic value of e^^ to be 1, and the factor of 

jj-^ numerator is there because of using the global symmetry to set e°"^ = 1. 
Let us assume that e^^ goes to 1 at infinity. Then the above simply states that what 
we have done is take some bare quantities ggeom and /Xgeom and combined them into 
the global symmetry invariant ggeomyUgeom = (e°'^|r=ooggcom) (e""'^|r=oo/^geom)- The 
second expression is made out of global symmetry invariants (e.g [e°'^\r=oo(lgeom)) , 
and so can define scales in the geometry. These are what we use to define the 
physical chemical potential /Xgeom, and why we will find e°"^\r=oo = 1 a convenient 
representative from all solutions related by the global symmetry. 

With the above considerations, and the perturbative results near the horizon, we 



may use the equations (26)-(28) and numerically integrate them from the horizon to 
a large value of r, where the functions settle to their AdS asymptotic values. These 
values at infinity furnish the required information to determine T and /x for the black 
brane. 



3 Results and Discussion 

Here we discuss the results found from numerically integrating the equations of mo- 
tion. In Figure [l| we give log-log plots of the dimensionless entropy density {'iGd+2s) 
versus the dimensionless temperature T for fixed fi = 1 and various dimensions. As 
seen in Figure [l| the entropy density is smooth and monotonic in T for a wide range 
of z and there are no possible discontinuous phase transitions associated with going 
from T <C /t to T ^ /t. The same behaviour is observed for the other dimensions 
2 < d < 9 not plotted in Figure [Tj Furthermore, we observe the correct asymptotic 
behaviour: the slope approaches d/z in the Lifshitz-like regime (T <^ /i); and the 
slope approaches d in the AdS regime {T ^ fi). 

We had thought from our earlier work |48] that the behaviour m. d = 2 might 
follow from the relatively low dimension of the dual model. However, the results 
of Figure [T] show that this is not the case. We see no evidence for a discontinuous 
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(a) d = 3. 



(b) d = 5. 
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Figure 1: The plots of In {AGd+2s) versus In(LT) for fixed ft = 1. Figures (la), (lb) 



(Ic) and (Id) correspond to d = 3, 5, 7, and 9, respectively. The different curves 
in each plot correspond to different values of a, with a = 4 green (solid), a = 2 
magenta (long-dashed), a = 1 cyan (dot-dashed), and a = 0.75 blue (dashed), a is 
related to z via a = \/2dl (z — 1). 
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phase transition or thermal instabihty regardless of the dimension. To be sure of 
these results required careful numerical examination of the dilaton at infinity to 
ensure that it was indeed smooth and approaching a constant value, as is appropriate 



for AdS. This is because the expression (25) for (j) has explicit r dependence, and 



so rounding errors will eventually occur for sufficiently large r. The powers of r 
involved are dimension dependent, so the rounding effect is more severe the higher 
the dimension. 

Next we turn to the question of measuring the number of degrees of freedom for 
the Lifshitz-like theories. We note first that the entropy density of the system in the 
Lifshitz-like limit is of the form 

AGa+2S = c{z,d){L^iY{^ ' , (47) 

where c(z, d) is a function of the dimension and the critical exponent. We have 
chosen /i = 1 in the plots shown in Figure [T| and so the above expression reduces to 
AGd+2S = c{z, d) T'^/^. Therefore, if we fit the ln(T) -C ln(/i) portion of the graphs in 
Figure [l] with a line, the intercept value of this line gives ln(c(2;, d)). This coefficient 
is a direct measure of the degrees of freedom of the theory. Figure [2] depicts the 
behaviour of ln(c(z, d)) versus ln(2;) for various dimensions. 

We can get a qualitative picture of how these must behave in the large z limit. 
First, note that as z becomes large, the slope of the T <^ fi part of the graph, 
d/z, goes to zero. The intercept, therefore, approaches a fixed value, and so c{z^d) 
must asymptote to a constant number. This is a trend we do observe in Figure [2j 
Therefore, c{z^ d) only depends on the dimension of the theory at large z. 

We cannot compare the values of c{z, d) we found for different z and different d, 
since they correspond to different theories. Nonetheless, it would be interesting to 
compare the values of c{z, d) we have found here for model S with the same quantities 
in different models. We leave this to future work. 

In the following, we will compare the large z behaviour of the graphs of Figure 
[2] with expectations from the dual field theory, using a simplified model. We start 
by considering a {2^T^Y volume box of some medium, and consider the excitations 
above the ground state. For such a case, the "legal" wave numbers may be written 



(48) 



v ■■ J 

Further we assume Lifshitz scaling symmetry and rotational symmetry in the spatial 
dimensions (broken only by the presence of the "box" ) . The only consistent relation 
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Figure 2: The plot of ln(c(z,(i)) as a function of \n{z) for fixed value of /t = 1. 
The different curves correspond to different dimension with d = 3 black (solid), 
d = 5 brown (long-dash), d = 7 red (dot-dash), and d = 9 coral (dash). Curves 
for d = 2,4,6,8 behave similarly. First of all, notice that the curves' intercept at 
ln(2;) = is given by the value In ((47r)'^/((i + l)"^) which is non-monotonic in d. 
Secondly, we find it interesting that the tails become flat out at 2; — > oo and that the 
large z behaviour for various d is monotonic in d. Given these two facts, it may be 
no surprise that the curves \n{c{z,d)) for fixed d are generically non-monotonic for 
small values of z, and that in fact they cross each other. 



between the energy of a state and the wave numbers above is 

= (v^)^/i^-^5 = i^'Yf^j^J , (49) 

where the powers of /i (an intensive parameter) appear by dimensional analysis. The 
coefficient 6 will in general be a function of z and d; for the time being we will leave 
this implicit, because we expect that 6 should remain order 1. Indeed, /i is supposed 
to be the relevant scale for new physics, and a large 6 would make this identification 
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suspect. 

Proceeding as in [48] , we assume that the occupation number of excitations above 
the ground state is e~^^^J^' (e"^"^"). This is merely a convenience since the function 
remains unspecified. The only content is that the occupation number is some func- 
tion of the ratio of the temperature to the energy of the state. An important physical 
assumption has gone in at this point. Essentially we have a bath of particles, with 
some number density n. We are assuming very small temperatures so that this den- 
sity may be considered very large, i.e. even though we are exciting a certain number 
of them above their ground state, the "same number" remain unexcited. Therefore, 
we are treating this number density of particles as an inexhaustible bath. This essen- 
tially means that the number of possible excitations is approximately infinite, and 
there is no effective chemical potential for adding new ones. The thermodynamics 
is governed by some collective mode, the "phonons", rather than the constituents of 
the material, the "molecules." 

Taking the above functional form of the occupation numbers, we may write 

E = -^F, F = 5^^(e-^-)- (50) 

n 

Calculating F, and approximating the sum as an integral, we find 



z {d/Sfiy/^ J h 



n 



-7 (e-^) (51) 



and define, therefore. 



z J J n 
Finally, we read that the energy is 



F(-)= [ n'/^^T{e-^). (52) 



E = ^^-A^T'^'^'f(^). (53) 



^2 (5^)rfA 

Note that the constant i'^ = Vd/{2TtY. Hence, we may take derivatives of the above 
by T at constant volume in a straightforward way. Therefore 
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Happily, this exhibits the T''/^ behaviour that we see in the T <^ fx hmit of the gravity 
system. Using the above expressions, one can show the expected thermodynamic 
relation 

e^j^Ts, (55) 

where we have divided by the volume to give the expressions in terms of the energy 
density S and entropy density s. The above energy that we have arrived at is the 
"dynamical" energy: it cannot be directly compared with the energy in AdS. The 
absolute energy in AdS will have a constant offset, due to the finite charge density and 
chemical potential of the background. Since this is some fixed "zero point" energy for 
low temperatures, it does not come into the above considerations in the effective field 
theorjj^ For a direct comparison to the numerical setup we would need to subtract 
off the asymptotic contribution: we do not do this here, and content ourselves with 
seeing the correct dependence of s on T for T <^ /i, which then implies the above 
expression up to an additive constant. 

Now let us recall that we have done the above for one type of particle. We must 
multiply this by the number of species to make a meaningful comparison to the 
gravity setup. In fact, the measure of the number of species is given by something 
like Ns ~ L"'/ (4Grf_|_2) for the gravity setup so the total entropy above gets multiplied 
by this amount. Finally, dividing by volume {27i£Y, we find that the total entropy 
density should be 



4G«. = 4G.„^=*(.,d)-/M.,<;)^^(i,.rF^;j [^f, (56) 

where on the right hand side the function k,{z, d) measures the number of species 
exactly: Ns = K{z,d)L'^ / {AGd+2)- 

One possible function of interest is a simple exponential suppression, J-'(x) = x, 
i.e. a Boltzmann distribution. This is in some sense a "dilute phonon" limit, where 
the excitation modes do not interact. In such a situation we find 

4G.,,. = i(.,d)-*'.(.,d)||i^±^(L^)^r(^) (I)'", (57) 

and the above is a good unitless measure of the entropy density. We can, of course, 
combine the two functions 5 and k together into one function that measures the 
number of degrees of freedom. 

^This is akin to ignoring the rest mass of a non-relativistic material when computing the energy. 
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We are now in a position to make a qualitative comparison to the graphs in Figure 
[2] First, we assume that k{z, d) and 6{z, d) are always order 1 numbers. Hence, in 
the large z hmit, these must asymptote to some fixed value, depending only on d. 
Therefore, all that is left to analyze is the behaviour of the factor F (^) {d + z)/ z^. 
In the large z limit F {d + z)/ z"^ — )■ ^ and so we see that the whole expression 
(57) goes to a 2; independent value, depending only on d. Therefore, the graphs in 



Figure |2j qualitatively agree with the large z behaviour of equation (57), given the 
restriction on 5{z^d) and K,{z,d). 



Acknowledgements 

This work was funded by a grant from NSERC of Canada. The work of IGZ was 
also supported by an Ontario Graduate Scholarship. 



References 

[1] J. M. Maldacena, Adv. Theor. Math. Phys. 2, 231 (1998) [Int. J. Theor. Phys. 
38, 1113 (1999)], |arxiv:h ep-th/9711200| S. S. Gubser , I. R . Klebanov and 
A. M. Polyakov, Phys. Lett. B 428, 105 (1998) |arxiv:hep-th/9802109|. E. Wit- 
ten, Adv. Theor. Math. Phys. 2, 253 (1998) | arxiv:hep-th/9802150j . 

[2] Y. Oz and J. Terning, Nucl. Phys. B 532, 163 (1998) [arxiv:hep-th/9803167|. 
S. Kachru and E. Silverstein, Phys. Rev. Lett. 80, 4855 (1998) |arxiv:hep-J 
fth/9802183]. A . E. La wrence, N. Nekrasov and C. Vafa, Nucl. Phys. B 533, 
199 (1998) |arxiv: hep-th/ 9803015|. I. R. Klebanov and E. Witten, Nucl. Phys. 
B 536, 199 (1998)lM^hep-th/9807080|. S. Benvenuti, S. Franco, A. Hanany, 
D. Martelh and J. Sparks, JHEP 0506, 064 (2005) |arxiv:hep-th/0411264|. 

A. Butti, D. Forcella and A. Zaffaroni, JHEP 0509, 018 (2005) |arxiv:hep^ 
|th70505220|. S. Franco, A. Hanany, D. Martelh, J. Sparks, D. Vegh and 

B. Wecht, JHEP 0601, 128 (2006) |arxiv:hep-th/050521 1]. 

[3] I. R. Klebanov and M. J. Strassler, JHEP 0008, 052 (2000) [arxiv:hep-J 
th/0007191|. J. M. Maldacena and C. Nunez, Phys. Rev. Lett. 86, 588 (2001) 
' [arxiv:hep-th/000 8001] . T. Sak ai and S. Sugimoto, Prog. Theor. Phys. 113, 843 
(2005) |arxiv:hep- th70412141|. T. Sakai and S. Sugimoto, Prog. Theor. Phys. 
114, 1083 (2005) |arxi v:hep-th /0507073|. 



18 



[4] S. W. Hawking and D. N. Page, Commun. Math. Phys. 87, 577 (1983). A. Cham- 
blin, R. Emparan, C. V. Johnson and R. C. Myers, Phys. Rev. D 60, 064018 
(1999) |arxiv:hep-th/9902170|. M. Cvetic and S. S. Gubser, JHEP 9904, 024 



(1999) [ arxiv:hep-th/9902195j . 
[5] D. T. Son, Phys. Rev. D 78, 046003 (2008) [arXiv: 0804. 3972] [hep-th]]. 



[6] K. Balasubramanian and J. McGreevy, Phys. Rev. Lett. 101, 061601 (2008) 
| arXiv:0804.4053 [hep-th]]. 

[7] T. Albash and C. V. Johnson, J. Phys. A 43, 345405 (2010) farXiv:0907.5406j 
[hep-th]]. P. Basu, J. He, A. M ukherjee and H. H. Shieh, Phys. Rev. D 82, 
044036 (2010) |arXiv:0908.1436l [ hep-th]]. S . S. Gubser, F. D. Rocha and P. Ta- 
lavera, JHEP 1010, 087 (2010) [arXiv:0911.3632 [hep-th]]. S^. Hartnoll and 
D. M. Hofman, Phys. Rev. B 81, 155125 (2010) |arXiv70912.00l)8 [cond-mat.str- 
el]]. S. A. Hartnoll, J. Polchinski, E. Silverstein and D. Tong, JHEP 1004, 120 
(2010) | arXiv:0912.106 1 [hep-th]]. T. Faulkner, N. Iqbal, H. Liu, J. McGreevy 
and D. Vegh, |arXiv~003.1728 [hep-th]. B. H. Lee, D. W. Pang and C. Park, 



JHEP 1007, 057 (2010) |arXiv:T0'06.1719 [hep-th]]. L. Y. Hung, D. P. Jatkar 



and A. Sinha, Class. Quant. Grav. 28, 015013 (2011) |arXiv:1006.3762 [hep- 
th]]. C. Hoyos-Badajoz, A. O'Bannon and J. M. S. Wu, JHEP 1009, 086 (2010) 
farXiv: 1007.0590 [hep-th]]. 



[8] A. Volovich and C. Wen, JHEP 0905, 087 (2009) |arXiv:0903.2455| [hep-th]]. 



[9] S. A. Hartnoll, Class. Quant. Grav. 26, 224002 (2009) |arXiv:0903.3246| [hep- 
th]]. 

[10] S. Sachdev, arXiv:1002.2947 [hep-th]. 



[11] C. P. Herzog, J. Phys. A 42, 343001 (2009) [arXiv: 0904. 1975 [hep-th]] 



[12] A. Ada ms, K. B alasubramanian and J. McGreevy, JHEP 0811, 059 (2008) 
[arXiv:0807.1111' [hep-th]]. C. P. Herzog, M. Rangamani and S. F. Ross, JHEP 
0811, 080 (2008) |arXiv:0807.1099 [hep-th]]. J. Maldacena, D. Martelh and 
Y. Tachikawa, JHEP 0810, 072 (2008) |arXiv:0807.1100 [hep-th]]. D. Yamada, 
Class. Quant. Grav. 26, 075006 (2009) |arXiv:0809.4928 [hep-th]]. L. Mazzu- 



cato, Y. Oz and S. Theisen, JHEP 09 04, 073 (2009) |arXiv:0810.3673 [hep 



th]]. A. Donos and J. P. Gauntlett, JHEP 0903, 138 (2009) [arXiv:0901.0818j 
[hep-th]]. N. Bobev, A. Kundu and K. Pilch, JHEP 0907, 107 (2009) 
|arXiv:0905.0673. [hep-th]]. A. Donos and J. P. Gauntlett, JHEP 0907, 042 



19 



(2009) tarXiv:0 905.109 8 [hep -th]]. A. Donos and J. P. Gauntlett, JHEP 0910, 
073 (2009) |arXiv:090 7!l76Tl [hep-th]]. E. Imeroni and A. Sinha, JHEP 0909, 
096 (2009) laTXiv: 0907.18921 [hep-th]]. 

[13] S. Schafer-Nameki, M. Yamazaki and K. Yoshida, JHEP 0905 (2009) 038 



|arXiv:0903.4245| [hep-th]]. 



K. Copsey and R. Mann, |arXiv:1011.3502| [hep-th]. 

J. Blaback , U. H. Danielsson and T. Van Riet, JHEP 1002, 095 (2010) 
tarXiv:100ll945| [hep-th]]. 



K. Balasubramanian, K. Narayan, JHEP 1008, 014 (2010) [arXiv: 1005.3291 
[hep-th]]. 

H. Singh, larXiv: 1009.0651] [hep-th]]. 



H. Singh, I arXiv:1011.6221| [hep-th]]. 

A. Donos and J. P. Gauntlett, JHEP 1012, 002 (2010) |arXiv:1 008.2"062l [hep- 
th]]. 

R. Gregory, S. L. Parameswaran, G. Tasinato and I. Zavala, JHEP 1012, 047 
(2010) [ arXiv: 1009.34451 [hep-th]]. 

A. Donos, J. P. Gauntlett, N. Kim and O. Varela, JHEP 1012, 003 (2010) 
(jtfXiv: 1009.38051 [hep-th]]. 

S. Kachru, X. Liu and M. Mulligan, Phys. Rev. D 78, 106005 (2008) 
farXiv: 0808. 1725] [hep-th]]. 



M. Taylor, arXiv:0812.0530| [hep-th]. 

A. Adams, A. Maloney, A. Sinha and S. E. Vazquez, JHEP 0903, 097 (2009) 
|arXiv:0812.0166 [hep-th]]. R. G. Cai, Y. Liu and Y. W. Sun, JHEP 0910, 080 
(2009) |arXiv:0909.2807 [hep-th]]. Y. S. Myung, Y. W. Kim and Y. J. Park, 
[arXiv:0910.4428 [hep-th]. M. C. N. Cheng, S. A. HartnoU and C. A. Keeler, 
JHEP 1003, 062 (2010) |arXiv:0912.2784ii [hep-th]]. E. Ayon-Beato, A. Garbarz, 



G. Giribet and M. Hassaine, JHEP 1004, 030 (2010) |arXiv: 1 001.2361 [hep-th]]. 
[25] C. M. Chen and D. W. Pang, JHEP 1006, 093 (2010) (aFXiv:1003.5064j hep-th]]. 



20 



[26] U. H. Danielsson and L. Thorlacius, JHEP 0903, 070 (2009) [arXiv:0812.5088l 
[hep-th]]. 



[27 
[28 

[29' 
[30 

[31 

[32 

[33 

[34 
[35 
[36 
[37; 
[38 

[39 

[40 

[41 
[42 



R. B. Mann, JHEP 0906, 075 (2009) |arXiv: 0905 .11361 [hep-th]]. 

G. Bertoldi, B. A. Burrington and A. Peet, Phys. Rev. D 80, 126003 (2009) 
tarXiv:0905.3183| [hep-th]]. 

D. W. Pang, |arXiv:0905.2678 | [hep-th]. 

G. Bertoldi, B. A . Burrington and A. W. Peet, Phys. Rev. D 80, 126004 (2009) 
|larXiv :0907.4755| [hep-th]]. 

M. H. Dehghani and R. B. Mann, Phys. Rev. D 82, 064019 (2010) 
I arXiv: 1006.3510 [hep-th]]. 

K. Balasubramanian and J. McGreevy, Phys. Rev. D 80, 104039 (2009) 
|arXiv:0909.026 3 [hep-th]]. 

E. J. Brynjolfsson, U. H. Danielsson, L. Thorlacius and T. Zingg, J. Phys. A 
43, 065401 (2010) |arXiv:0908.2611 [hep-th]]. M. H. Dehghani and R. B. Mann, 
JHEP 1007, 019 (2010) |arXiv: 1004.4397 [hep-th]]. 



S. F. Ross and O. Saremi, JHEP 0909, 009 (2009) tarXiv:0907.1846| [hep-th]]. 



Y. Liu and Y. W. Sun, JHEP 1007, 099 (2010) |arXiv: 1006.2726 [hep-th]]. 



K. B. Fadafan, arXiv:0912.4873 [hep-th]. 



S. J. Sin, S. S. Xu and Y. Zhou, |arXiv:09 09.4857' [hep-th]. 

M. Ca doni, G. D'AppoUonio and P. Pani, JHEP 1003, 100 (2010) 
(a^iv:0912.3520| [hep-th]]. 



F. Apri le, S. Franco, D. Rodriguez- Gomez and J. G. Russo, JHEP 1005, 102 
(2010) farXiv: 1003.44871 [hep-th]]. 

C. Charmousis, B. Gouteraux, B. S. Kim, E. Kiritsis and R. Meyer, JHEP 1011, 
151 (2010) I arXiv: 1005.4690 [hep-th]]. 



B. H. Lee, S. Nam, D. W. Pang and C. Park, arXiv: 1006.0779] [hep-th]. 



E. Perlmutter, ^arXiv: 1006.2124 [hep-th]. 



21 



[43] S. A. Hartnoll, P. K. Kovtun, M. Muller and S. Sachdev, Phys. Rev. B 76, 
144502 (2007) |arXiv: 0706. 32 15| [cond-mat .str-el] ] . 

[44] 



[45 

[46; 
[47 



[49 
[50 



W. Y. Wen, [arXiv: 1009.3952 [hep-th]. 



V. G. M. Puletti, S. Nowling, L. Thorlacius and T. Zingg, |arXiv:1011.6261 
[hep-th] . 



S. A. Hartnoll and P. Petrov, arXiv: 1011.6469 ] [hep-th]. 

K. Goldstein, S. K achru, S. Prakash and S. P. Trivedi, JHEP 1008, 078 (2010) 
[larXiv: 091 1.35861 [hep-th]]. 

G. Bertoldi, B. A. Burrington and A. W. Peet, Phys. Rev. D82, 106013 (2010) 
arXiv: 1007. 1464 [hep-th]. 

S. W. Hawking and G. T. Horowitz, Class. Quant. Grav. 13, 1487 (1996) 
|arxiv:gr-qc/9501014j . 

A. Batrachenko, J. T. Liu, R. McNees, W. A. Sabra and W. Y. Wen, JHEP 
0505, 034 (2005) [ arxiv:hep-th/0408205] . 



22 



